Color superconductivity with determinant interaction in strange quark matter 
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<N ! 

We investigate the effect of six fermion determinant interaction on color superconductivity as well 
as on chiral symmetry breaking. Coupled mass gap equations and the superconducting gap equation 
are derived through the minimisation of the thermodynamic potential. The effect of nonzero quark 
- antiquark condensates on the superconducting gap is derived. This becomes particularly relevant 
for the case of 2-flavor superconducting matter with unpaired strange quarks in the diquark channel. 
While the effect of six fermion interaction leads to an enhancement of u-d superconductivity, due to 
nonvanishing strange quark-antiquark condensates, such an enhancement will be absent at higher 
densities for u-s or d-s superconductivity due to early (almost) vanishing of light quark antiquark 
condensates. 
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I. INTRODUCTION 



Cold and dense quark matter can be a color superconductor. At sufficiently high baryon densities, when nucleons get 
converted to quark matter, the resulting quark matter is expected to be in one kind or the other of the many different 
Oh, possible color superconducting phases at low enough temperatures 0. The rich phase structure is essentially due to 
the fact that the quark-quark interaction is not only strong and attractive in many channels, but also many degrees 
of freedom become possible for quarks like color, flavor and spin. Studying the properties of color superconducting 
phases in heavy ion collision experiments seems unlikely in the present accelerators as one cannot avoid producing 
large entropy per baryon in heavy ion collisions and hence cannot produce the dense and cold environment that is 
needed to support the formation of superconducting phases. However, in the future accelerator facility planned at GSI 
for compressed baryonic matter experiments, one possibly can hope for observing fluctuations signifying precursory 
phenomena of color superconducting phase [2| . 

On the other hand, it is natural to expect some color superconducting phases to exist in the core of compact stars 
where the densities are above nuclear matter densities and temperatures are of the order of tens of KeV. However, to 
consider quark matter for neutron stars, color and electrical charge neutrality conditions need to be imposed for the 
bulk quark matter. Such an attempt has been made in Ref.Q as well as in Ref. 4], where the lighter up and down 
quarks form the two flavor color superconducting (2SC) matter while the strange quarks do not participate in pairing. 
A model independent analysis was done in Ref. Q that is valid in the limit m s << /i and A ~ m} s j[i, where, A is 
the pairing gap and fi is the quark chemical potential. It has been shown, based upon the comparison of free energies 
that such a two flavor color superconducting phase would be absent in the core of neutron stars Q . Within Nambu 
Jona-Lasinio (NJL) model in Ref. |J] it has been argued that such conclusions are consistent except for a small window 
in density range where superconducting phase is possible. There have also been studies to include the possibility of 
mixed phases p) of superconducting matter demanding neutral matter on the average. Later, it was observed that 
the imposition of neutrality conditions lead to pairing of quarks with different Fermi momenta giving rise to gapless 
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modes |al3- Within a Nambu Jona Lasinio model, the two flavor superconducting quark matter (2SC) was shown to 
exhibit gapless modes (g2SC) arising due to the difference in the Fermi momenta of the pairing quarks, when charge 
and color neutrality conditions are imposed. Superconducting quark matter with unpaired strange quarks (2SC+s) 
was shown to exhibit gapless superconductivity (g2SC+s) within a window of about 80 MeV in baryon chemical 
potential 0. Temperature effects on the gapless modes were also studied for the two flavor quark matter in Ref. 0,0. 

A variational ap pro ach was used to study the chiral symmetry breaking as well as color superconductivity for the 
quark matter liAllOj . The calculations were carried out for the N JL model with the minimisation of the free energy 
density to study which condensate will exist at what density. Charge neutrality conditions were introduced through 
the introduction of appropriate chemical potentials. We note here that the possibility of diquark condensates along 
with quark-antiquark condensates has also been considered in Ref.s 0, 0, 0, Q] . 

In the present work, we investigate the effect of six fermion determinant interaction on color superconductivity in 
Nambu JonaLasinio model. Since it is not possible to have a model with same symmetries as QCD with four quark 
operators alone, the so called t'Hooft term is added to the usual four fermion interaction which is of the form 

Cdet ~ det f $(l + j 5 )ip] + det f $(l - j 5 )ip] (1) 

where, determinant is in the flavor (u,d,s) space. This term respects chiral symmetry and breaks U(l) axial symmetry 
as in QCD. 

While st udying the dense quark matter, the effect of such a term in quark antiquark channel has been widely 
used 0, 0, Hal ITq . However, for color superconductivity such a term is o nly considered at the four fermion level. 
The effect of six fermion interaction has been considered recently in Ref. [17]]. This, however considers a different 
effective six fermion interaction than as given in Eq.Q). We shall treat both chiral symmetry breaking as well as color 
superconductivity in the same footing using the determinant interaction as in Eq.Q. With the variational ansatz 
considered, we show that such a term gives a nonzero contribution to the free energy and at higher densities, it is of 
similar or even larger magnitude than the contributions from the quark-antiquark condensates. 

We organize the paper as follows. In the following subsection, we discuss the ansatz state with the quark-antiquark 
as well as the diquark condensates 0] . In section [HJ we consider the Nambu Jona-Lasinio model Hamiltonian with 
the determinant interaction and calculate its expectation value with respect to the ansatz state to compute the energy 
density as well as the thermodynamic potential. We minimise the thermodynamic potential to calculate all the ansatz 
functions and the resulting mass as well as superconducting gap equations here. In section HTT1 we discuss the results 
of the present investigation. Finally we summarise and conclude in section llVl 



A. The ansatz for the ground state 

We shall use the notations and conventions of Ref. 0, El and recapitulate briefly the construction of the variational 
ansatz for the ground state. We take it as a squeezed coherent state involving quark-antiquark as well as diquark 
condensates as given by U ITo| 

|0) =U d \vac) =U d U Q \0). (2) 

Here, Uq and Ud are unitary operators creating quark-antiquark and diquark pairs respectively. Explicitly, the 
operator, Uq is given as 

U Q = cxp (J q 0i (k)i((T ■ k)hi{k)q 0l {k)dk - h.c^j . (3) 

In the above, q° (q°) is the two component particle annihilation (antiparticle creation) opertor of the four componet 
quark field operator ip. The superscript '0' refers to the fact that these operators correspond to "free" Dirac fields 
including, in general, a current quark mass. The quark-antiquark condensate function hi(k) is a real function of 
|k| which describes vacuum realignment for chiral symmetry breaking for quarks of a given flavor i. We shall take 
the condensate function h(k) to be the same for u and d quarks and h^ik) as the chiral condensate function for the 
s-quark. Clearly, a nontrivial hi(k) shall break the chiral symmetry. Summation over three colors and three flavors 
is understood in the exponent of Uq in Eq. ^ . Similarly, the unitary operator Ud describing diquark condensates is 
given as 



U d = exp(B\ - B d ) 



(4) 
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where, B\ is the pair creation operator as given by 



l (k)V/(k)gi b r (-k)t e „- e3ab + gT(k)r/i(k)gi b r (-k)e ye3 aJ dk 



(5) 



In the above, i,j are flavor indices, a, b are the color indices and r(= ±1/2) is the spin index. The operators <z(k) are 
related to g°(k) through the transformation q(k) = Uqq (\l)Uq 1 . As noted earlier we shall have the quarks of colors 
red and green (a=l,2) and flavors u,d («=1,2) taking part in diquark condensation. The blue quarks (a=3) do not 
take part in diquark condensation. Note that we have assumed the 'condensate functions' /(k) to be independent 
of flavor color indices. We give a post-facto justification for this to be that the function depends upon the average 
energy and average chemical potential of the quarks that condense and is the same for red up and green down or 
green up and red down quarks, when color isospin is unbroken. 

Finally, to include the effects of temperature and density we next write down the state at finite temperature and 
density |fi(/3,/i)) taking a thermal Bogoliubov transformation over the state |Q) using thermofield dynamics (TFD) 
as described in Rcf.s 18, 191. We then have, 



where Up^ is given as 



with 



B t (/3,/x) = J [^(k) t 0-(k,/3, A1 )g;(k)t+g;(k)0 + (k,/?, A1 )|;(k)' 



dk. 



(6) 



(7) 



(8) 



In Eq.JSJl the ansatz functions 0±(k, ft, fx) will be related to quark and antiquark distributions, and, the underlined 
operators are the operators in the extended Hilbert space associated with thermal doubling in TFD method. In Eq.© 
we have suppressed the color and flavor indices on the quarks as well as the functions 0(k, f3, fi). 

All the functions in the ansatz in Eq.© are to be obtained by minimising the thermodynamic potential. We shall 
carry out this minimisation in the next section. 



II. EVALUATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 

As discussed earlier we shall consider here 3-flavor Nambu Jona Lasinio model including the determinant interaction 
given by equation ft|. The Hamiltonian is then given as 



H = ij)\-ia • V + 7 °m)V - G s J2 [OM^) 2 - (ip^ 5 X A ip) 2 ] 



A=0 



G D (^ 7 5 eeV C ')(^ C 7 5 «V) + K [dei/0(l + l5 )^\ + det f $(l - 75 



(9) 



where ip l,a denotes a quark field with color 'a', (a = r,g,b) and flavor l V, (i = u, d, s) indices. The matrix of 
current quark masses is given by m=diag/(?n u , m^, m s ) in the flavor space. As noted earlier, we shall assume isospin 
symmetry with m u —md- In Eq.Q, X A , A — l,---8 denote the Gellman matrices acting in the flavor space and 

A = 1/ as the unit matrix in the flavor space. The four point interaction term ~ G s is symmetric in 

SUv(3) x SUa(3) x LV(1) x ^a(I)- In contrast, the determinant term ~ K which for the case of three flavors 
generates a six point interaction which breaks ?Za(1) symmetry. If the diquark and the mass terms are neglected, the 
overall symmetry is SUv(3) x SUa(3) x C/y(l). This spontaneously breaks to SUy (3) x J7y(l) implying the conservation 
of the baryon number and the flavor number. The current quark mass term introduces additional explicit breaking 
of chiral symmetry and the axial flavor current is not completely conserved. 

The third term in Eq.@ describes a scalar diquark interaction in the color antitriplet and flavor antitriplet channel. 
We have not included, in the present investigation, a pseudoscalar term with the same coupling which leads to 
Goldstone mode condensation. Such a form of Lagrangian can arise e.g. by Fiertz transformation of a four point 
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current current interaction having quantum numbers of single-gluon exchange |20| . In that case the diquark coupling 
Gd is related to the scalar coupling as Gd = 0.75G S . 

Using the variational ansatz state in Eq.© one can calculate the expectation values of various operators 0. We 
evaluate the expectation values 



and, 



where, 



(n(ft M )l#?(k)$WW.M)> = ^A^k^/z^k-k'), 
(0(^, M )|^(k)-t^( k ')|0(^, M )) = S ij <5 ai A!^ (k, 0,/j,)6(k k'), 



A£(k,/3, M ) = i [l± (Fi a (k) - F ia (k)) ± ( 7 °cos^(k)) + a ■ ksin^(k))(l - F ia (k) - F?(k)) 



(10) 



(11) 



(12) 



In the above, V>(k) is the Fourier transform of @. The effect of the diquark condensates and their temperature 
and/or density dependences are encoded in the functions F m (k) and F-£ a (k) given as 



F™(k) = sin 2 0!°(k) + sin 2 / cos (T jb (k) (l - 5 



ra3\ 



(13) 



and, 



ff(k) = sin 2 6»; a (k) + sin 2 f x cos 29f jb (k) (l - 6 a3 ) . (14) 

1 — sin 2 0± — sin 2 6± with i ^ j and a ^ b. The <5 a3 term indicates that the third 
color does not take part in diquark condensation. Further, we have introduced the notation for the quark-antiquark 



condensate functions as a 'shift' from their vaccum values as (j>i(k) = </>?(k) — 2/ij(k), with cot </>°(fe) = mi/ (k 2 +mf). 
We also have for diquark operators, 



(fi(/3,M)|C a (xX(0)|O(/3, M )) 
(^(/3, M )|C at (x)V4 fct (0)|^(/3,M)> 



1 



(2tt)3 

1 

w 



e lk ' x n a 4i 6 (k,/3,M)rfk, 



(15) 



where, 



n° ji, (k,/3, M ) 



ij ~3ab 



e' J £ 



7° cos 



S ia ' jb (k) cos 



cc • k sin 



7 5 C, 



(16) 



and, 



pi ojb (k,/3,/i) 



S ia ' jb (k) cos 



+ 7 cos 



a • k sin 



(17) 



Here, C = i^ 2 j° is the charge conjugation matrix (we use the notation of Bjorken and Drell) and the functions S(k) 
and A(k) are given as, 



and, 



S^ ab (k) = sin2/(k) cos26» wjb (k,/?,^) +sin2/i(k) cos 2(9! h ajb (k, /?, /i), 
A ij '' o6 (k) = sm2/(k)cos26>! ajh (k,/3,M) - sin2/i(k) cos 20™ jb (k, 0, /x), 



(18) 
(19) 
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These expressions are used to calculate thermal expectation value of the Hamiltonian and compute the thermodynamic 
potential. 

Let us first concentrate on the contribution of the determinant term to the energy expectation value. When 
expanded the determinant term will have six terms, each involving three pairs of quarks of different flavors. These 
are to be 'contracted' among themselvs in all possible manner, while taking expectation values. Further, while 
considering quark-antiquark condensates, it is clear that the contractions of the same color will be dominant over 
that with different colors by a factor N c , where iV c is the number of colors. With the present case of two flavor 
superconductivity, this leads to the fact that out of the six terms in the expansion of the determinant, only two 
terms that are proportional to the strange quark-antiquark condensate (ss) will be dominant over the rest. These 
latter four terms are suppressed at least by a factor N c . The dominant two terms are the ones involving contraction 
of strange quark-antiquarks having the same color. This simplification arises because we are considering only u-d 
superconductivity here. For the rest of our calculations, we take contributions of these two terms only. Explicitly 

these two terms are given as ~ Sa(*O a s) (uO a u)(dO a d) — (uO a d)(dO a u) , where we have written ,with a = ±, 

(5 ± = (1 ± 75). Using, Eq.s (|15|l -Ea. I|17|). it is starightforward to show that, both the terms in the square bracket 
give identical contribution to the expectation value when 'contracted' diquarkwise, except, that the contribution of the 
second term is of opposite sign as compared to the first term due to the flavor antisymmetric nature of the expectation 
values. The determinant term contribution is then given as 

V det = +K{det f ${l+^}+det f [${l-^ 5 )i)]) 

= -8KlP(ll + 2lPlP). (20) 

Here, 

3 

^° = ~^v < >=E / j^™^-***-*?) ( 21 ) 

is proportional to the quark-antiquark condensate for i-th flavor, and, 

I D = i(Tp Q 7 5e y 'e 3Q V b ) 

= — i-3 f dkcos( ^ 1 ~ ^ ) [sin2/(k)(l - sin 2 0i - sin 2 6> 2 ) + sin2/ 1 (k)(l - sin 2 0\ - sin 2 6\)] (22) 

is the diquark condensate. It is then straightforward to calculate the expectation value of the Hamiltonian. This can 
be written as 

e = (H) = T + V s + V D + V de t, (23) 

the various terms arising from the kinetic, the scalar, the diquark and the determinant interaction terms of the 
Hamiltonian respectively. Explicitly, the kinetic energy part in Eq. (JHJ is given as 

T = (0(AM)IV' + Ha-V + 7V-)^|^(AM)) 
2 3 3 



,••) -, 3 EE / dk d k l sin ^ + m * cos -F ia -Ff), (24) 

^ 71 > i=l a=l 

where, F la and F{ a are given by Ea.s fT3*jl and Ea. lfHjl . 

The contribution from the scalar interaction term in Eq. Q turns out to be 

V S = -G s (fi(/3,M)I E [«^) 2 " (h 5 ^) 2 } = -8G S £ /f, (25) 

A=0 i=1.3 



where, I\ is given in Eq.JSlJ. 
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Similarly, the contribution from the diquark interaction from Eq.(|5Jl to the energy density is given as 

V D = -Gz,(r!(/3,/i)|(^ 7 5 e b V C )W C 7 5 «V)^(/3,M)> = -16G D I% (26) 

with Id defined in Ea. (|22|l and the contribution from the determinant term is as given in Ea. (|20() . 

To calculate the thermodynamic potential we shall have to specify the chemical potentials relevant for the system. 
Here we shall be interested in the form of quark matter that might be present in compact stars older than few minutes 
so that chemical equilibration under weak interaction is there. The relevant chemical potentials in this case are the 
baryon chemical potential /i.b(= 3/i g ), the chemical potential [Le associated with electromagnetic charge, and the two 
color electrostatic chemical potentials /Z3 and /ig. The chemical potential is a matrix that is diagonal in color and 
flavor space, and is given by 

fHj,ab = i^ij + QijfJ-E)$ab + (Q3abfJ-3 + QsdbPs)&ij • ( 27 ) 

Demanding the color superconducting ground state to be invariant under SU(2) C gauge group, we can choose the 
chemical potential /13 to be zero. 

The total thermodynamic potential, including the contribution from the electrons, is then given by 

n = T + Vs + V D -{»N)-~s + Sl e (28) 

where, we have denoted 

(jtN) = (^fx ij<ab ^ b ) = 2£> to #. (29) 

i,a 

In the above, zx M is the chemical potential for the quark of flavor i and color a, which can be expressed in terms of 
the chemical potentials zx, [le and zxs using Ea. Q27|l . Further 



7 *° = (2^ / dk{Fm ~ Fn 



(30) 



is proportional to the number density of quarks of a given color and flavor. The thermodynamic potential for electrons 
is given as 



/4 ( 2 t 2 

12^ 2 V Ml 



«e = -^Z5 ( 1 + 2 ^7t)- (3 1 ) 



Here, the electron mass is assumed to be zero which suffices for the system we are considering. 
Finally, for the entropy density for the quarks we have |18| 



j^s Yl J dk { sin2 ° l - ln sin2 6 - + cos2 6 - ln cos2 0i ~ 



( 27 

+ sin 2 6™ ln sin 2 0™ + cos 2 6™ ln cos 2 6^ a ) . (32) 

Now the functional minimisation of the thermodynamic potential £1 with respect to the chiral condensate function 
hi(k) leads to 

. ... m i + 8G.ll i) +8K\e» k \Is^I^ +4KI D S i3 M l 
cot ^(k) = j^j = |^|-. (33) 

The last term in the numerator indicates the explicit dependence of the strange quark-antiquark condensate function 
on the light diquark condensate Id- 

Substituting this back in Ea. l21|) yields the mass gap equation as 



Mi =m, + ^ / _L £ (1 _ F io - F?)dk + 8K\e ijk \I^I^ + 4KI D 6* 3 , 

£ * a =l,3 



(34) 
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with, ej = \J\s? + Mf being the energy of the constituent quarks of i-th flavor. Note that for the two flavor super- 
conductivity as considered here, the strange quark mass is affected explicitly by the superconducting gap given by 
the last term on the right hand side Ea. l|34[l . Of course, there is implicit dependence on the gap in the second term 
through the functions F and F\ (given in Eq.s Ijl3|l and l|14|l). Further, when chiral symmetry is restored for the 
light quarks i.e., when the scalar condensate for the nonstrange quarks vanishes, still, the determinant term gives 
rise to a density dependent dynamical strange quark mass. Such a mass generation is very different from the typical 
mechanism of quark mass generation through quark-antiquark condensates |l7j . 

Next, the minimisation of the thermodynamic potential Q with respect to the diquark and di-antiquark condensate 
functions /(k) and /i(k) yields 



and 



tan2/(k) 



tan 2/i (k) 



4(2G D + Klf ] ) 
l-p 



A ,<k 

; COS( 



4(2G D +Kli 3) I D ) A 



e + p 



P 



■ cos(- 



(35) 



(36) 



where, we have defined the superconducting gap A = 4(2Gd — KI^)Id, with ljj as defined in Ea. (|22|l . Further, 
in the above e = (ei + £2)/2, p — (p U r + Mdg)/2 = (p U g + Pdr)/2 = A 1 + Ms/6 + /is/vo. It is thus seen that the 
diquark condensate functions depend upon the average energy and the average chemical potential of the quarks that 
condense. We also note here that the diquark condensate functions depends upon the masses of the two quarks which 
condense through the function cos ((<f>i — </>2)/2). The function cos0i = Mj/ej, can be different for u,d quarks, when 
the charge neutrality condition is imposed. Such a normalisation factor is always there when the condensing fermions 
have different masses as has been noted in Ref. 0] m the context of CFL phase. 

Substituting the solutions for the condensate functions given in Eq.s (|33J) . Ij35(l and (|36|l in the expression for In in 
Eq. (|22|l . we have the superconducting gap equation given by 



(2tt) 3 J \ 2 J 



1 



cos 26L 



cos 29 



(37) 



In the above, Q± = A 2 cos 2 ((0i - fa)/2) + ^± = e ± p, ,and, cos26»± = 1 - sin 2 6%. - sin 2 6»^. 

Finally, the minimisation of the thermodynamic potential with respect to the thermal functions 0±(k) gives 

sin 2 6^ - 1 



exp(/3w5?) + 1 



(38) 



Various uj la, s (i, a = flavor, color) are explicitly given as 



UJ 
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.12 



gj± + 5 e ± 6u 



0J4 



(39) 



uj 2 } = uf = u>± - S e ± = u>i (40) 
for the quarks participating in condensation, and, 

c4° = e l ±fi ia (41) 

for the blue quarks which do not participate in superconductivity. Here S e (= (ei — £2)/2) is half the energy difference 
of the two quarks which condense, and, <5 M = (pu — M22V2 = /i_e/2, is half the difference of the chemical potentials 
of the two condensing quarks. Further u>± have been already defined after Ea.Q37|). Note that when the charge 
neutrality conditions are not imposed, all the four quasi particles taking part in diquark condensation will have the 
same energy u)_. It is clear from the dispersion relations given in Ea. (|41(l that it is possible to have zero modes, i.e., 
u! la = depending upon the values of S e and 8^. So, although we shall have nonzero order parameter A, there will 
be fcrmionic zero modes or the gapless superconducting phase |2jJ, |22J . 
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Now using these dispersion relations, the mass gap equation Eq. H34|) and the superconducting gap equation Eq. <|37|) . 
the thermodynamic potential (Eq.lJSSJ) becomes 



a, 



(42) 



In the above, Q q is the contribution from the quarks and is given as 



a, 



(2tt) 3 
4 



dk 



•\/k 2 + to 2 (w_ + uj + ) 



+ 



(3(2 
4 



[ — Jdk [log(l + exp(-/3w!?) + log(l + exp(-/30] 



(2tt) 3 
2 



m—u,d 

dk 



\Jk 2 + m 2 - -(ex + e 2 ) 



/3(2 
6 



(2tt) 3 
2 



l —J2[dk [log(l + exp(-/3( ei - M l3 )) + log(l + exp(-/3(e, + ^ i3 ))] 
dk \^Jk? + m 2 - yjk? + M 2 



Jdk [log(l + exp(-/3(e 3 - p 3a )) + log(l + exp(-/3(e 3 + M 3a ))] 



^ 2 -) 3 a tT,3 

16GD/I, + 8G S Y l( s l)Z + 32X1^1^1^ + IQKlfh 



2 

D I 



(43) 



i=l,3 



where, are given in equations (j39|l and (|40[1 . The contribution of the electron 57 e to the total thermodynamic 
potential is already given in Ea. lp?T|l . The first two lines in Ea. lf^ correspond to the contributions from the quarks 
taking part in the condensation while the third and fourth lines correspond to the contribution from the two light 
quarks with the blue color. The next two lines are the contributions essentially from the strange quarks. The last line 
corresponds to the terms which mix up the flavors and also the diquark and quark-antiquark condensates. In fact, 
the last two terms in the expression of thermodynamic potential are the contributions arising due to the determinant 
interaction. 

Thus the thermodynamic potential is a function of four parameters: the three mass gaps and a superconducting gap. 
These are calculated through minimisation of the thermodynamic potential, subjected to the conditions of electrical 
and color charge neutrality. The electric and color charge neutrality constraints are given as 



and, 



Qi 



Qs 



l 



l 



2p 



0. 



0, 



(44) 



(45) 



respectively. In the above p la = (ip la ^ip' la ) = 2I l v a (i, a not summed) and T l v a is as given in Ea. H3U|) . Further, 
p 1 = J2 a =i 3 P ia ■ The thermodynamic potentials (Eqs. (|42|) and (|43|I % 1. the charge neutrality conditions (Eqs. (|44|l and 
(|45in . the mass gap equation (Eq.||!2J)) an d the superconducting gap equation (Eq.||23)) constitute the basis of the 
numerical calculations that we shall discuss in the next section. 



III. RESULTS AND DISCUSSIONS 



For numerical calculations we have taken the values of the parameters of the NJL model as follows. The coupling 
constants G s , Go have the dimensions of [Mass] -2 while the six fermion coupling has a dimension [Mass] -5 . To 
regularise the divergent integrals we use a sharp cut-off, A in 3-momentum space. Thus we have all together six 
parameters, namely the current quark masses for the nonstrange and strange quarks, m q and m s , the three couplings 
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Fig. 1-a 
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Fig.l-b 



FIG. 1: Gap parameters when charge neutrality conditions are not imposed. Fig. 1-a shows the gaps at zero temperature as a 
function of quark chemical potential. Fig. 1-b shows the same when the determinant interaction is not taken into account. Both 
the figures correspond to nonzero value for the current quark masses m u =5.5 MeV—md- Solid curve refers to masses of u and 
d quarks, the dotted curve refers to mass of strange quarks and the dashed curve corresponds to the superconducting gap. 



G s , Gd, K and the three- momentum cutoff A. For simplicity, we shall also take Gd to be 0.75G S , as may be expected 
from Fierzing a current current interaction. For the rest of the parameters we choose A = 0.6023 GeV, G S A 2 = 1.835, 
KA 5 = 12.36, m q = 5.5 MeV and ra s = 0.1407 GeV as has been used in Ref. 23]. After choosing m q = 5.5 MeV, 
the remaining four parameters are fixed by fitting to the pion decay constant and the masses of pion, kaon and vj '. 
With this set of parameters the mass of r\ is underestimated by about six percent. With this parametrization, the 
constituent masses of the light quarks turn out to be Mi = 0.368 GeV for u-d quarks, while the same for strange 
quark comes out as M s — 0.549 GeV, at zero temperature and zero density. 

It is, however, relevant here to comment regarding the choice of the parameters. There have been different sets of 
parameters by other groups also [2^, l26l [27J for the three flavor NJL model. Although the same principle as above is 
used e.g. in Ref j2^|, the resulting parameter sets are not identical - in particular, the dimensionless coupling KK° 
differs by as large as 30 percent as compared to the value used here. This discrepancy lies on different treatment of the 
?y' meson. Since NJL model docs not confine, and because of the large mass of the rf meson (m ?) /=958 MeV), it lies 
above the threshold for qq decay with an unphysical imaginary part for the corresponding polarization diagram. This 
is an unavoidable feature of NJL model and leaves uncertainty which is reflected in the difference in the parameter 
sets by different groups. Within this limitation regarding the parameters of the model, however, we proceed with the 
above parameter set which has already been used in the study of the phase diagram of dense matter in Ref. fl5| as 
well as in the context of equation of state for neutron star matter in Ref. p4j . 

Let us begin with the discussions of results when the charge neutrality conditions are not imposed. At zero 
temperature, the behaviour of the gap parameters as functions of quark chemical potential are displayed in Fig.l- 
a. We may point out that these solutions for the gaps correspond to the solutions for which the thermodynamic 
potential is minimised. In fact, in general, for certain values of the chemical potential, there can be several solutions 
of the gap equations particularly near the critical chemical potential. We have chosen the ones which minimise the 
thermodynamic potential. 

At low chemical potentials fi q < /ii ~ 350 MeV, the diquark gap vanishes and the masses of the quarks stay at their 
vacuum values. The entire region below \i q = [i\ corresponds to vacuum solution and has zero baryon number. At 
fi q = H\ a first order phase transition takes place and the system is a two flavor color superconductor. The diquark 
gap jumps from zero to about 95 MeV at this point. At the same point, the masses of u and d quarks drop from 
their vacuum values of 370 MeV to 50 MeV. The baryon number density also jumps from zero to 0.42 fm -3 at this 
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chemical potential. Because of the flavor mixing six Fermi interaction, this cross over transition for the light quarks 
is reflected also as dropping of the strange quark mass at this chemical potential from its vacuum value of 549 MeV 
to about 470 MeV. Although this transition is not a first order transition due to the nonzero current quark masses, 
there exist metastable phases. The masses of the light quarks in these metastable phases are the nontrivial solutions 
of the mass gap equations, but, have higher thermodynamic potential as compared to the solutions corresponding to 
stable phases which are shown in Fig.l-a. With the increase in fi q , the superconducting gap increases until it reaches 
a maximum at fi q ~ 475 MeV. Beyond this point, as the quark chemical potential is increased, the mass of strange 
quark drops again. Due to the determinant interaction this leads to a drop in superconducting gap as may be obvious 
from Eq. I|37l) . For comparison we have plotted in Fig. 1-b the masses and the superconducting gap for the case where 
the determinant coupling is zero, while all other couplings remain unchanged. In this case, it may be noted that 
the (nearly) vanishing of light quark masses does not have any effect on the strange quark mass. The strange quark 
mass starts to drop only when quark chemical potential is larger than the vacuum value of the mass of the strange 
quark as only then the finite density contributions become nonzero in the mass gap equation for the strange quark. 
In contrast, with nonzero coupling K, because of the flavor mixing, the dropping of strange quark mass starts for 
quark chemical potentials smaller than the vacuum mass of the strange quark as shown in Fig. 1-a. Since the charge 
neutrality condition depends very sensitively on the strange quark mass, determinant coupling therefore can have 
important consequences in deciding the phase structure. 

We would also like to note that, to evaluate the thermodynamic potential, we have to solve self consistently the 
three mass gap equations fEa.l|34|l) and the superconducting gap equation fEa. (|37[l 'l. These equations are all coupled. 
However, some simplification occurs for higher densities if we take the current quark masses to be zero. This leads to 
solving only two coupled gap equations - the strange quark mass gap equation and the superconducting gap equation 
for densities when chiral symmetry is restored for the light u,d quarks. Although the nature of chiral symmetry 
transition changes from a cross over to a first order transition, the values of the constituent masses of the quarks 
do not change very much. Henceforth we shall limit our discussions to this case when current quark masses for the 
light quarks are taken to be zero. In Fig. 2 we show the results of such a calculation for zero temperature when the 
charge neutrality conditions are not imposed. The general behaviour of the graphs are similar to those in Fig. 1. The 
important difference is that, the chiral transition is a sharp first order transition at zero temperature. The constituent 
quark masses at zero densities for up (down) and strange quarks are 354 MeV and 546 MeV as compared to 368 MeV 
and 549 MeV respectively, when current quark masses of u,d quarks are taken to be nonzero. The critical chemical 
potential turns out to be 335 MeV for the case of nonzero K. We also note that as the quark chemical potential 
increases beyond 480 MeV, the superconducting gap starts decreasing as the strange quark-antiquark condensate 
decreases so that the effective diquark coupling also decreases, as may be clear from Eq. l|3TJ) . 

There have been calculations where the effect of the determinant interaction is retained for mass gap calculations, 
but the effect of such an interaction is not included in the superconducting gap equation 0, 0, We plot in 

Fig. 3 the superconducting gap including the effect of determinant interaction for the u-d superconductivity gap 
for the case when the charge neutrality conditions are not imposed. As may be seen, including the effect of strange 
quark-antiquark condensate through the determinant interaction leads to an enhancement of the superconducting 
gap and this can be as large as a twenty five percent for chemical potentials of about 400 MeV. This is because the 
'effective' coupling for superconductivity increases from Go to Gd — K(ss)/4. At higher densities however, this effect 
diminishes as the maginitude of the quark-antiquark condensate itself decreases. Such effects can play an important 
role in deciding which phase can occur at what density as we come down in the density. In particular, such kind of 
enhancement will be there for u-d superconductivity, but will be absent for u-s or d-s superconducting gaps since the 
light quark-antiquark condensate vanishes much earlier. 

We next extend our discussion to the case when the charge neutrality conditions (Qe=0=Qs) are imposed. We 
compute the thermodynamic potential numerically as follows. For given values of the quark chemical potential, /i q , 
and the electric and color charge potentials, he and /is, the coupled mass gap equations and the superconducting gap 
equations are solved self consistently. The values of /i£ and /is are varied so that the charge neutrality conditions 
fEa. (|44|l and Eq.|@SJ) are satisfied. The resulting solutions are then used in Ea. l|43|l to compute the thermodynamic 
potential. In doing so, we also check for existence of multiple solutions of the gap equation and if they exist, we choose 
the solution which has the least value for the thermodynamic potential. In Fig^J we show the dependence of the 
superconducting gap on quark chemical potential when charge neutrality condition is imposed. The superconducting 
gap starts becoming nonzero for [i q greater than 350 MeV. At this point the chiral symmetry for the light quarks is 
restored. This has also its effect on the strange quark mass which drops from its vacuum value of about 546 MeV 
to 470 MeV. The superconducting gap increases smoothly with /j, q until /i g attains a value of around 430 MeV. At 
this point the gap jumps from 80 MeV to 106 MeV and then increases slowly upto a maximum value of about 107.5 
MeV at n q = 450 MeV. Beyond this value of fx q , the magnitude of the strange quark-antiquark condensate decreases 
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FIG. 2: Same as in Fig. 1, but with zero current quark masses 




FIG. 3: Superconducting gap as a function of quark chemical potential. Solid line corresponds to including the effect of strange 
quark-antiquark condensate. The dashed line corresponds to the case when this effect is not included. 



leading to a drop of superconducting gap through the determinant interaction. 

As the quark chemical potential increases, strange quarks help in maintaining the charge neutrality conditions. We 
might observe here that the strange quark mass starts decreasing already at \i q = 425 MeV, when charge neutrality 
conditions are imposed. This may be compared with [i q ~ 480 MeV, when charge neutrality conditions are not imposed 
(see Fig. 2a). As we have already mentioned, there could be multiple solutions of the gap equation and whichever 
has the least free energy is the stable solution. It could so happen that a solution of the mass gap equation which 
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FIG. 4: Superconducting gap as a function of quark chemical potential for the charge neural quark matter 

is metastable when charge neutrality conditions are not imposed, can be the only solution when charge neutrality 
conditions are imposed. 

The pairing between the quarks with charge neutralty conditions correspond to stressed pairing i.e. the pairing 
of the quarks of different species which differ in their Fermi momenta |2^| . This gives rise to possibility of a gapless 
superconducting phase, the QCD analogue of Sarma phase |29| . In the present case, between the quark chemical 
potentials fi q ~ 350 MeV and n q ~ 425 MeV, the system is in the gapless phase. The number densities of u and d 
quarks participating in the condensation are plotted in Fig. Vs the quark chemical potential is increased beyond 
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FIG. 5: Number densities of u quarks (solid) and d quarks (dot-dashed) participating in superconducting phase. 
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425 MeV the solution for the gap jumps to a higher value of about 106 MeV almost similar to the case when charge 
neutrality conditions are not imposed. This corresponds to the usual BCS solution. In this case, the number densities 
of u quarks and d quarks participating in the condensation are the same. The charge neutrality conditions however 
are maintained by the blue colored u,d quarks, the electron as well as the strange quarks. One essential effect of 
including strange quarks is that the electron density starts decreasing for chemical potentials greater than the strange 
quark mass and the strange quarks carrying the negative charge maintain electric charge neutrality condition. This 
has the effect that the lowest excitation energy u>2 — ui + 8^ in the BCS pairing case becomes large due to both the 
large value of the gap as well as due to the small magnitude of the electron chemical potential. 

It may be worthwhile to mention here that the gapless modes in superconductivity were known theoretically in 
the context of superconducting matter with finite momentum |2l) as well as in condensed matter systems with 
magnetic impurities j2^. Recently it has been investigated for cold fermionic atoms [2H |^|. Gapless modes in 
the context of quark matter has been first proposed in Ref. for color flavor locked matter. However, this 

corresponded to a metastable phase. Gapless modes in neutral quark matter were first emphasized for the two flavor 
color superconductivity in Rcf. @ and for the 2SC+S quark matter in Ref. 7J. Stable gapless modes for color flavor 
locked matter were first proposed in Ref. and have been confirmed in Ref. |l5j| in a more general structure for the 
gap. The temperature dependence of the CFL gapless modes has been studied in Ref.(3j|. We might mention here 
that the effects of nonzero strange quark mass and the charge neutrality condition on superconducting gaps have also 
been studied within a Ginzberg Landau approach in Ref. |34| . 



IV. SUMMARY 



We have analysed here the effect of flavor mixing t'Hooft six fermion interaction both for chiral symmetry breaking 
as well as two flavor superconductivity in NJL model. The mehod is a variational one with an explicit construct for the 
trial state including quark-antiquark as well as diquark condensates. The determinant interaction affects explicitly 
both the u-d superconducting gap as well as the mass gap for the strange quarks. For the strange quarks, a density 
dependent mass arises even when the quark-antiquark condensates vanish. Such type of generation of dynamical mass 
is entirely distinct from the typical mechanism of spontaneous chiral symmetry breaking through quark-antiquark 
condensates 0. Further, the u-d superconducting gap gets enhanced as inclusion of such a term effectively increases 
the diquark coupling by an amount proportional to the strange quark-antiquark condensate. 

We have focussed our attention here to the two flavor superconducting phase, with unpaired strange quarks. The 
variational method adopted can be directly generalised to include color flavor locked phase and one can then make 
a free energy comparison regarding the possibility of which phase would be thermodynamically favorable at what 
density. In case of CFL phase, such an enhancement of superconducting u-s (d-s) condensates will not be there as 
the corresponding (dd)((uu)) will vanish much earlier than (ss) condensates, when the density is increased. Since 
strange quark mass is a sensitive parameter in maintaining the charge neutrality condition, this will be important 
while comparing free energies of different phases of charge neutral quark matter. 
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